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I ports energy from the left edge of a linear ion crystal to the right edge 

OO ' by local operations and classical communication at a speed consider- 

. ably greater than the speed of a phonon in the crystal. A probe qubit 

I is strongly coupled with phonon fluctuation in the ground state for a 

' short time, and it is projectively measured in order to obtain informa- 

■ tion about this phonon fluctuation. During the measurement process, 

phonons are excited by the time-dependent measurement interaction, 
K/i \ and the energy of the excited phonons must be infused from outside 

j_j ' the system. The obtained information is transferred to the right edge 

^ of the crystal through a classical channel. Even though the phonons 

excited at the left edge do not arrive at the right edge at the same 
time as when the information arrives at the right edge, we are able to 
soon extract energy from the ions at the right edge by using the trans- 
ferred information. Because the intermediate ions of the crystal are 
not excited during the execution of the protocol, energy is transmitted 
in the energy transfer channel without heat generation. 



1 Introduction 



A trapped ion system is expected to be a promising candidate for devel- 
oping quantum computers [1]. Recently, the experimental development of 
trapped ions has undergone a significant technical advancement [2]. In ad- 
dition, this technology has been applied to quantum teleportation |3] and 
quantum error correction [1]. The successful application of trapped ion sys- 
tems to the abovementioned tasks proves that these systems have a great 
potential for use in other fascinating quantum tasks, which have not yet 
been experimentally executed. In this study, we analyse quantum energy 
teleportation with cold ions in a linear trap, that is, a linear ion crystal. Re- 
cently, quantum energy teleportation (QET) has been proposed this 
protocol is able to transport energy simply by local operations and classical 
communication. Energy can be effectively teleported without breaking any 
physical laws including causality and local energy conservation by measuring 
zero-point oscillation at one site in the entangled ground state of a many- 
body quantum system and providing the measurement result to distant sites. 
The key point is that there exists quantum correlation between local fluc- 
tuations of different sites in the ground state. Therefore, the measurement 
result of local fluctuation in some site includes information about fluctuation 
in other sites. By selecting and performing a suitable local operation based 
on the transferred information, the zero-point oscillation of a site away from 
the measurement site can be more suppressed than that of the ground state, 
yielding negative energy density. Here, the origin of energy density is fixed 
such that the expectation value vanishes for the ground state. It is well known 
that such regions with negative energy density are created by superposing 
energy eigenstates quantum mechanically. Even if a system comprises a re- 
gion with negative energy density, other regions with positive energy density 
still exist, and therefore, the total energy of the system remains non-negative. 
Discussions on negative energy density of non-relativistic systems are given 
in ^ and those for negative energy density of relativistic systems are given 
in [H]. While performing the above local operation to create negative energy 
density in the system, the surplus energy is transferred from quantum fluc- 
tuations to external systems and can be harnessed. Though the protocols of 
energy teleportation can be implemented for other systems such as quantum 
fields [5] and spin chains [6l[71[8], linear ion crystals can be effectively used for 
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QET. Due to the Coulomb interaction and the harmonic potential for linear 
trapping, the ions interact with each other so strongly that entanglement for 
energy teleportation is easily generated in the phonon ground state. It is also 
possible to POVM measure the local fluctuation of phonons in the linear ion 
crystal by coupling the phonon modes with the internal energy levels of each 
ion via a laser field [2]. 

In this study, we use the linear trap models that are extensively analysed 
by James [10] and consider a protocol of QET. In what follows, N cold ions, 
which are strongly bound in the y and z directions but weakly bound in a 
harmonic potential in the x direction, form a linear ion crystal called the 
QET channel. The first ion that stays at the left edge of the crystal is the 
gateway of the QET channel where energy is input. The A^-th ion that stays 
at the right edge of the crystal is the exit of the QET channel where the 
teleported energy is output. We select two suitable internal energy levels of 
the gateway ion and regard them as the energy levels of a probe qubit in 
order to measure the local phonon fluctuation. The probe qubit is strongly 
coupled with the phonon fluctuation in the ground state for a short period 
of time via a laser field, and it is projectively measured in order to obtain 1- 
bit information about phonon fluctuation. During the measurement process 
at the gateway, phonons are excited in the system by the time-dependent 
measurement interaction, and hence, the energy of the excited phonons must 
be infused from outside the system. In the measurement models used in this 
study, the kinetic energy of the gateway ion increases after the measurement; 
however, the kinetic energy of other ions and the potential energy of all the 
ions remain unchanged. This infusion of energy is regarded as energy input 
at the gateway of the QET channel. The obtained information is transferred 
through a classical channel from the gateway point to the exit point. In 
principle, the speed at which this information is transferred is equal to the 
speed of light, which is considerably greater than that of phonon propagation 
in the ion crystal. It is emphasized that even when the phonons excited at 
the QET gateway do not arrive at the exit point, the information still arrives 
at the exit point. Because the measurement process at the gateway is local, 
the state of the exit ion before the phonons arrive is locally the same as 
the ground state. Surprisingly, however, we are able to soon extract energy 
from the exit ion by using the transferred information. Performing a local 
operation that depends on the measurement result suppresses the zero-point 
oscillation of the phonon fluctuation at the exit ion, yielding negative energy 
density around the exit. The surplus energy of the fluctuation is transferred 
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from the phonon system to the external systems including the device, in 
order to perform this operation. Here, it should also be emphasized that 
without the measurement result, we cannot extract energy from the exit ion 
in the local ground state by performing an arbitrary local quantum operation 
on the ion. Because the intermediate ions of the crystal are not excited 
during the protocol execution, this QET ensures energy transportation in 
the channel without the generation of heat. The non-negativity of the total 
Hamiltonian of phonons ensures that the amount of the teleported energy is 
less than the amount of the input energy. The amount of energy required for 
transferring the information is negligibly small in principle as compared to 
the amount of the teleported energy, because the information can be simply 
transported by low-energy physical carriers such as electromagnetic waves 
with long wavelengths. 

The paper is organized as follows. In section 2, the linear ion crystal 
model developed by James [10] is reviewed. In section 3, QET for the model 
is analysed. In the last section, the summary and discussion are given. In 
this study, we adopt a unit h = 1. 



2 Linear Ion Crystal Model 



Following James [10] , let us consider N ions with charge Ze and mass m that 
are strongly bound in the y and z directions but weakly bound in a harmonic 
potential in the x direction. The position of the n-th ion is denoted by 
where the ions are numbered from left to right. The Hamiltonian is expressed 
as 

^ m ^ 1 ^ Z'^f? 

£t 2 ^ £t 2 ^ 2 ^^^^ \xn{t) - ^ ^ 

where the dot denotes time derivative, and u is the trap frequency, which 
characterizes the strength of the trapping harmonic potential in the x direc- 
tion. In the case of cold ions, we can approximate the position of the n-th 
ion by 

where Xn^ is the position of the ion in the ground state, and g„(t) is a 
small displacement describing the phonon modes in the ion crystal. The 
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equilibrium positions in the ground state are determined by the following 
equation derived from Eq. ([T]): 

4^ ( (0) (o)\^ ^ / (0) (0)^2 ■ ^ ' 

Let us introduce a scale length given by 

/ = 

and rescale the position variables as 



mz/2 



Ur. 



I 

Then, Eq. ([2]) is rewritten in a dimensionless form as 

n-l ^ ^1 

First, it is pointed out [lOj that Eq. ([3]) can be analytically solved for the 
cases with = 2 , 3. When N = 2 , the solution of the above equation can 
be expressed as 

2 2 
1\5 



When = 3 , the solution can be expressed as 

For greater than 3, m„ were numerically solved by James. The Hamiltonian 
of phonons in the crystal is derived by expanding Eq. ([1]) in terms of qn{t) 
and by considering only bilinear terms; it is expressed as follows: 

^1 ^1 
^ = Yl ^mu^Ann'qnqn' - Eg, (4) 

n=l n,n'=l 
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where the zero-point energy Eg is subtracted from the original form in order 
to make the lowest eigenvalue of H zero, and a real symmetric matrix A^n' 
is defined by 



^ 1 

A„„ = 1 + 2 ^ ^ (5) 

for n = n' and 



n"=l 



A , 



1 



for n 7^ n'. The eigenvectors 6n (A; = 1, 2, ■ ■ ■ , A^) of Ann' determined by 



N 



n'=l 



are real. In addition, the eigenvalues /i^ of Ann' are real and non-negative 
|10j . The eigenvectors are numbered in the order of increasing eigenvalues 
and normalized as 



TV 



n=l 



The first and second eigenvectors are analytically obtained. For k = 1, the 
eigenvector is given by 



'N 

with its eigenvalue /ii = 1. For k = 2 , the eigenvector is given by 



6?) 



Yln'=l ^n' 



with its eigenvalue /i2 = 3. Higher eigenvalues and eigenvectors are numeri- 
cally solved and listed in table 2 of [TU]. Next, let us introduce the normal 
modes of phonons as 



N 



n=l 



Then, H is diagonalized for the modes as 

N 



/I 1 \ 

k=l ^ ^ 



where Pk is the conjugate momentum of Qk and Pk{t) = mQ]^{t). Let us 
introduce phonon creation and annihilation operators as 



, 1 D . 



These operators satisfy the following commutation relation: 

^kk'- 



The ground state \g) of the phonon system is defined by 

ak\g) = 0, (6) 

and simultaneously satisfies 

H\g) = (7) 

due to the subtraction of the zero-point energy Eg. From Eq. ([7]), it is 
verified that if is a non-negative operator, i.e. H > 0. By using the normal 
modes, it is easy to solve the Heisenberg operator for the displacement of the 
n-th ion from the equilibrium position as 

N 



k=l 
N 

= E /^A^ (a,e-v^* - at e-^/^*) . (8) 



Its corresponding conjugate momentum operator is solved as 



N 



k=l 
N 



^^(,) (a.e-v^* + 4e-v^*) . (9) 

k=i 
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These Heisenberg operators can be computed from their Schrodinger opera- 
tors Qni^ ?n(0)) and Pn(=Pn(0)) as follows: 



Quit) = Wiiithw + X: W^S(t)Pn', 

n'=l n'=l 

N N 

n'=l n'=l 

where W!j[^,{t) (r — 1, 2, 3) are real symmetric matrices defined by 

N 

k=l 



fc=l 
N 



k=l 

The operators ?„(i) and describe the quantum motion of phonons in the 
crystal. Under real experimental conditions, the typical order of the largest 
frequency v^/Jin of the phonon oscillation is O(IO^) Hz for = 2 ~ 10, and 
the typical order of the ion crystal size is O(10~^) m. Therefore, the ratio 
of the phonon velocity to the light velocity is estimated to be 0(10"^) <C 1. 
This estimation ensures that the system can be treated non-relativistically. 

In the later discussion, it is convenient to introduce phonon coherent 
states. Let us consider a unitary operator given by 



U (a, P) = exp 



N 



n=l 



where a = {ai, ■ ■ ■ , ctAr) and /? = {Pi, ■ ■ ■ , P^) are real vectors. By performing 
U (a, P) on the ground state, a coherent state is obtained as follows: 



\{a,P)) = U{a,P)\g). 
This state is an eigenstatc of the operator expressed as 

1 



(10) 



Ofcl {a,p)) 



[a,P)), 
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where 



N 



Ak = J^^n^^n, 

n=l 
N 

Bk = J2bl^^Pn. 



n=l 



Solving Eqs. ([8]) and IQ with t = 0, \ can be exphcitly written as 



exp 



X exp 



1 ^ 



k'=l 



1 



■ N 



1 



Ay -i\ mv./Ji^Bk 



Ak Bk < 



\9). 



11) 



In addition, the inner product between two coherent states is calculated as 

' . N 

2 

(Ak - A[) - iJmu^kiBk - B[) . (12) 



exp 



X exp 



. N 

I 

1 ^ 



k=l 



^mv^k 



These formulas are used in the next section. 



Quantum Energy Teleportation 



In this section, we analyse a QET protocol by treating the linear ion 
crystal discussed in the previous section as a QET channel. The first ion 
at X = is the gateway of the QET channel where energy is input. The 
iV-th ion at a; = xf^ is the exit of the QET channel where the teleported 
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energy is output. First, let us perform a local POVM measurement of the 
phonon fluctuation in the ground state \g) at the gateway. The probe qubit 
used to measure the fluctuation is composed of two internal energy levels of 
the gateway ion in the same way as [2j. The measurement interaction in our 
model is given by 

Hm = 9it)ayG,, (13) 

where g(t) is a time-dependent real coupling constant, ay is the y component 
of the Pauli matrix, and Gi is a Hermitian local operator on the first ion 
defined by 

Gi = (P + \qi. (14) 

Here, and A are time-independent real coupling constants. Let us assume 
that the coupling constant g{t) does not vanish only during a very short 
period of time via the laser pulse field and that it can be approximated as 

g{t) = 6{t). 

The initial state of the probe qubit is assumed to be the up eigenstate of the 
z component of the Pauli matrix, az, given by 



1+) 



1 




After the switch off of the measurement interaction in Eq. (ITSll . we pro- 
jectively measure of the probe qubit. This POVM measurement can be 
described by two measurement operators [TT] given by 



M± = (d=|exp [-iayGi] |+), 
where |— ) is the down eigenstate of given by 



(15) 




1 



The operators M± are computed explicitly as 



M4 
M_ 



cosGi, 
sinGi. 



(16) 
(17) 
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In addition, they satisfy the following relations: 

5^MtM, = l, (18) 

s=± 

sMjMs = cos (2Gi) . (19) 

s=± 

The average state of the phonon system after the measurement is given by 

PM = J2M,\g){g\Ml 

s=± 

Before the measurement, the phonon system is in the ground state \g) and 
has no energy such that 

{9\H\9) = 0. 

After the measurement, the average energy of the phonon system is given by 

= Tr [Hpm] = Y,{g\MlHMs\g). (20) 

s=± 

This can be easily evaluated using Eq. f|T8l) . Taking into account that Mg 
commutes with the kinetic energy of the ions for n = 2 ~ and the potential 
energy of all the ions, it can be proven that only the kinetic energy p\/ (2m) 
of the gateway ion changes during the measurement. In order to evaluate 
this change, the following relation can be used: 

[pi, M±] = ±i\M^ (21) 

From Eq. (!2T!) . we can obtain the following relation given by 

MiplM± = mIm± {pI + A^) T 2iXM+M_pi. (22) 

Using Eqs. (fT5|) and (1221) . Ei^ is computed as 

^'^ = ^i^\^^ni;MM-{,\^\,}^^. (23) 

s=± 

This energy of the excited phonons must be infused from outside the system 
and is regarded as the energy inputted at the gateway of the QET channel. 
The obtained information s is transferred through a classical channel from 
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the gateway point at x = x\ to the exit point at x = x)^ . In principle, the 
speed at which the information is transferred is equal to the speed of light, 
which is considerably greater than that of phonon propagation in the ion 
crystal. It is emphasized that the phonons excited at the QET gateway do 
not arrive at the exit point when the information arrives at the exit point. 
Because the measurement process at the gateway is local, the state of the 
exit ion before the arrival of the phonons is locally the same as the ground 
state. Interestingly, however, we are able to soon extract energy from the 
exit ion by using the transferred s. Let us perform a local unitary operation 
dependent on s given by 



where ^ is a real parameter fixed later. It is possible to suppress the zero- 
point oscillation of the phonon fluctuation at the exit ion by performing the 
above operation. Neglecting time evolution of the system during the high- 
speed transfer of s, the average state of the system after the operation is 
given by 




(24) 




(25) 



s=± 



Using Eq. (120!) and a formula given by 




N 



n=l 



the average energy after the operation can be evaluated as follows: 



Ef = Tr [Hpf] = Y,^g\MlUlHUsM,\g) 



s=± 
N 




= E,^ - 9mu^ Yl ^Nn{g\qn '^^sMs 1 \g) + O'^^A^n- (26) 



n=l \s=± 



Here, we have used the fact that [Mj, g„] =0. Substituting Eq. (|T9|) into 
Eq. ^ yields 

Ef = ~97^ + e'^i, (27) 



11 



where rj and ^ are real coefficients given by 



N 



7] = mu' 



^ANniglQn cos {2Gi) \g), 



n=l 



^ = Ann- 



(28) 
(29) 



The expression of rj in Eq. fl25]) can be simphfied further as follows. First, it 
is pointed out that the following relation holds: 



cos(2Gi) \g) = -[e^^G^ + e-''''^] \g) 



-2i(j> —2i\q 



(30) 



It is observed that the two states | ± 2A) = e^^'^'^'^'^lg) in the above equation 
are the following phonon coherent states in Eq. ffTOj) : 

|±2A) = |((±2A,0,---0),(0,0,---0))). 

Therefore, the states | ± 2A) are the eigenstates of such that 



afcl ±2A) 



2Xb\ 



(fc) 



±2A). 



(31) 



By introducing a matrix A„„/ given by 



^ 1 

An.'=i^i:.^,(o)=$:^«) 

tit 



and using Eq. ([6]), Eq. ^ with t = 0, and Eq. ( !3TI) . the following relation 
is proved. 

(0|g„| ±2A) = ±z — (0| ±2A)Ai„. (32) 



mv 



In addition, Eq. ( fT2l) yields the following relation: 

" A2 



± 2A) = exp 



-A 



mi> 



11 



(33) 
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From Eq. (1501) . Eq. (15^ . and Eq. (15^ . the following relation is obtained: 



(fif|g„cos (2Gi) \g) 



X sin (20) 



mz/ 



exp 



mu 



-An 



A 



In- 



Substituting this relation into Eq. ( l28i) yields the final expression of as 
follows:. 

An 



-Az/ sin (20) exp 



mz/ 



TV 



n=l 



From this result, we can show that the coefficient t] does not vanish as long 
as the factor sin (20) does not vanish. This observation indicates that the 
operation in Eq. (l24l) ensures that Ep is smaller than Ein. In fact, Ep in 
Eq. ( 1271) can be minimized by considering the parameter 6 such that 
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Then, Eq. ( l27j) can be rewritten as 



2f 



E, 



Eir). 



(34) 



Because A^]^ in Eq. ([5]) is positive, ,^ is also positive. This implies that 
during the operation Us-, a positive amount of energy given by 



E. 



out 



4e 



A^ sin^ (20) 



2m A 



exp 



NN 



2X' 



mu 



An 



N 



nN 



n=l 



(35) 



is transferred from the phonon system to the external systems including the 
device system executing the operation Ug. For example, Eout for = 2 is 
analytically evaluated as 



E. 



Vs X' 



out 



4 2m 



sin^ (20) exp 



1 



1 + ^ — 
V3/ mu. 



X' 



This energy extraction is regarded as outputting the teleported energy from 
the exit of the QET channel. Figure 1 shows the QET protocol. Eq. (135!) 

can be expressed in terms of the input energy E^ ^= ^ j as follows: 



Eout = iNEin exp -C 



Eir^ 



sin^ (20) 



(36) 
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where the coefficients 7iv and (n are defined as 



In 
Cn 



Ann 
4An. 



n 

E 

n=l 



In^nN 



(37) 
(38) 



The maximum value of the output energy Eout with respect to (p is obtained 
by setting </ 



= ±j. From Eq. fl36l) . it is observed that i^o^i exponentially 
decays when the input energy Ein is considerably greater than the typical 
energy u of one phonon in the ion crystal. Therefore, Ein should be of the 
same order as this is ensured by choosing A = O (i/mz/) in Eq. ( |T^ . 
Further, it is observed that Eout rapidly decays when becomes large. In 
Figures 2 and 3, numerical results of In 77V and (n are plotted as functions 
of N. Though (n does not have a drastic iV-dependence, 'Jn decays approx- 
imately exponentially (oc e~^'^^). Hence, the output energy Eout behaves as 
^-i.inq j|- emphasized that without using the transferred s, we cannot 
extract energy from the exit ion on an average by performing an arbitrary 
local quantum operation on the ion. This is because the state of the exit ion 
is locally the same as the ground state. Due to the passivity property of the 
ground state, any local quantum operation independent of s on the exit ion 
gives energy to the crystal, generating an excited state, or it ensures that the 
state is unchanged with no energy gain. The non-negativity of H ensures 
that Ep > in Eq. (IMI) . Thus, Eout is not greater than Em- 



4 Summary and Discussion 



In this study, we have analysed a QET protocol for the linear ion crystal 
model developed by James [10]. At the gateway of the QET channel (at 
X = x^^^), the phonon fiuctuation in the ground state is POVM measured. 
The measurement operators are given by Eqs. 0161) and 0171) . During the 
measurement process, the input energy Ein in Eq- is infused from out- 
side the system. The kinetic energy of the gateway ion increases after the 
measurement; however, the kinetic energy of the other ions and the poten- 
tial energy of all the ions remain unchanged. The obtained information s 
by the POVM measurement is transferred through a classical channel to the 
exit point at x = x^^\ In principle, the speed at which the information is 
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transferred is equal to the speed of light, which is considerably greater than 
that of the phonon propagation in the ion crystal. Even when the phonons 
excited at the gateway point do not arrive at the exit point, we are able to 
soon extract the energy Eout in Eq- fl35|) from the exit ion by performing Us 
in Eq. fl2^ . Because the intermediate ions of the crystal are not excited dur- 
ing all the operations of the protocol, this QET involves the transportation 
of energy without the generation of heat in the channel. 

Thus far, even though the QET mechanism is an interesting phenomenon, 
it has not yet been experimentally verified. Clearly, the experimental ver- 
ification of QET requires fast performance of local operations and classical 
communication. Therefore, it may be better to perform these operations 
collectively using one physical carrier of information. To realize such a situ- 
ation, polarization of a laser pulse in a fibre as a probe qubit can be carried 
out instead of using the internal energy levels of the ion as the qubit. The 
laser pulse is coupled sequentially with the first and iV-th ions in the crys- 
tal during its propagation in the fibre. This situation is depicted in Figure 
4. To make the discussion more concrete, let us consider the creation and 
annihilation bosonic operators ^l/J(C) and ^'s(C) for one photon of the laser 
field with polarization s = ± in the fibre parametrized by a coordinate (. 
The fibre is connected between the initial point {( = (i)and the final point 
{( = (f) via the first ion of the crystal at C = Ci ^ind the A^-th ion of the 
crystal at C = Ca^; where Ci < Ci < Ca^ < C/- The operators \I/J(C) and ^^^(C) 
satisfy the following commutation relations: 



Let us assume that the initial state of the laser field is a pulse- wave coherent 
state with polarization s = + given by oc exp (^f fi{()'^\^{()d(^ |0), 

where /i(C) is a function with a support localized around ( = Q. Let us 
consider a free Hamiltonian of the fibre photon as 



5ss'6 (c - a , 



0. 





oo 



— oo 
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where c is the hght velocity and ^I/(C) is given by 



^+(0 
*-(C) 



The free evolution of the photon field is given by 



e^*^*^l(C)e-**^* = - ct). 



In order to simulate the protocol discussed in section 3, this laser field couples 
with the first and the iV-th ions by interactions given by 



where d is the length of the interaction region. The total Hamiltonian of the 
composite system is expressed as 



and it is independent of time. Hence, Htot is conserved in time. The initial 
state of the composite system is given by \g)g ® In this model, the 

evolution of the laser pulse induces effective switching of the interactions 
for the QET operations on the phonon system. The interactions in Eqs. 
f l39|) and fHOl) are active only when the pulses exist at their corresponding 
interaction regions. Based on the interaction in Eq. fl5^ . the information s 
about the phonon fluctuation is imprinted into the polarization of the laser 
pulse. The interaction in Eq. (HOj) gives a controlled operation gate of the 
exit ion by the value s of cr^ of the laser pulse. The energy of the laser pulse 
changes when the pulse passes through the ion regions. The initial energy 
of the pulse is denoted by Ei. In principle, the photon energy Ei is chosen 
to be independent of Eout- In order to perform a sensitive experiment for 
the detection of the QET effect, it is better to consider that the order of Ei 
is the same as that of Eout{= 0{'^nv))- After the interaction with the first 
ion in Eq. fl39l) . the pulse energy is decreased to E2 by the exciting phonons 
with energy Ei^ = Ei — E2 in the crystal. After the interaction with the 
A^-th ion in Eq. fHU]) . the pulse energy is increased to E3 by extracting the 




(39) 



(40) 



tot — 



H + Hm + Hlo + He 
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energy from the ion as the QET effect. The amount of the extracted energy 
is E,j.i,i = E3 — E2. Experimental observation of this change in the pulse 
energy leads to the verification of the QET mechanism. Detailed analysis of 
this model will be reported elsewhere. 
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Figure Captions 

Figure 1: Schematic diagram of QET for the hnear ion crystal. The 
trapped cold ions are denoted by circles. The POVM measurement defined 
by Eqs. f|T6|) and f|T7|) is performed on the ion at the left edge of the crystal 
to obtain the information s about phonon fluctuation in the ground state 
\g). During the measurement, energy Em is infused into the crystal. The 
measurement result s is transferred to the right edge of the crystal through a 
classical channel. The energy Eout is extracted by performing Us in Eq. ( 12^ 
on the ion at the right edge. 

Figure 2: In 77V is plotted as a function of N. 

Figure 3: (n is plotted function of N. 

Figure 4: Schematic diagram of the verification experiment of QET. The 
black line denotes a photonic fibre for a laser pulse that controls the switching 
of the interactions for the QET protocol. The fibre is parametrized by a 
coordinate ( and connected between the initial point {( = Q) and the final 
point {( = Q) via the first ion of the crystal at C = Ci ^ind the A^-th ion of 
the crystal at C = Ca^- The initial energy of the laser pulse is denoted by 
El. After the interaction with the first ion in Eq. fl39|) . the pulse energy is 
decreased to E2 by the exciting phonons with energy Ein = Ei — E2 in the 
crystal. After the interaction with the A^-th ion in Eq. fHUl) . the pulse energy 
is increased to E3 by extracting the energy from the ion as the QET effect. 
The amount of energy extracted is given by Eout = E3 — E2. 
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